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ABSTRACT 

The main aim of microlensing experiments is to evaluate the mean mass of massive 
compact halo objects (MACHOs) and the mass fraction of the Galactic halo made 
by this type of dark matter. Statistical analysis shows that by considering a Dirac- 
Delta mass function (MF) for the MACHOs, their mean mass is about that of a white 
dwarf star. This result is, however, in discrepancy with other observations such as 
those of non-observed expected white dwarfs in the Galactic halo which give rise to 
metal abundance, polluting the interstellar medium by their evolution. Here we use the 
hypothesis of the spatially varying MF of MACHOs, proposed by Kerins and Evans 
to interpret microlensing events. In this model, massive lenses with a lower population 
contribute to the microlensing events more frequently than do dominant brown dwarfs. 
This effect causes the mean mass of the observed lenses to be larger than the mean 
mass of all the lenses. A likelihood analysis is performed to find the best parameters of 
the spatially varying MF that are compatible with the duration distribution of Large 
Magellanic Cloud microlensing candidates of the MACHO experiment. 

Key words: galaxies: halos-dark matter. 



1 INTRODUCTION 

The rotation curves of spiral galaxies (including Milky Way) 
show that this type of g alaxies have dark halo structure 
iBorriello fc Salucci 200 ll) . One of the candidates for the 
dark matter in the Galactic halo may be massive compact 
halo objects (MACHOs). Paczyriski (1986) proposed a gravi- 
tational microlensing technique as an indirect way of detect- 
ing MACHOs. Since his proposal, many groups began moni- 
toring millions of stars of the Milky Way in the directions of 
the spiral arms, the Galactic bulge and the Large and Small 
Magellanic Clouds (LMC fc SMC) and detected hundreds 
of microlensin g candidates (Ansari 2 0041 : iDerue et al. 20011 
ISumi et al. 2003 : Afonso et al. 20 03). Looking in the direc- 
tion of the LMC and SMC (which are the most impor- 
tant for estimating Galactic halo MAHCOs), Experience 
de Recherche d'Objets Sombres (EROS)^ and MACHO^ 
observers found only a dozen of micro lensing candidates 
jLasseree et al. 200ol : lAlcock et al. 2000h . The interpreta- 
tion of LMC and SMC events is based on the statistical 
analysis of the distribution of the duration of the events. The 
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result of this analysis attribute a mean mass to MACHOs 
and their mass contribution in the Galactic halo. With the 
standard halo model, the mean mass of lenses is evaluated 
to be about half of that of the solar mass with a 20 per cent 
contribution in the Galactic halo mass. 
The results obtained by the analysis of LMC microlens- 
ing events, however, do not agree with other observations 
JCates fc Gvuk 200 ll) . Studying the kinematics of white 
dwarfs that have been discovered dOppenheimer et al. 200lh 
has shown that halo white dwarfs corresponds to 1 — 2 
per ce nt of the halo mass . Recent re- analysis of the same 
data i Spaena et al. 20041 : iTorres et al. 20021) shows that 
this fraction is an order of magnitude smaller than the value 
derived in Oppenheimer et al. (2001). On the other hand, 
if there were as many white dwarfs in the halo as suggested 
by the microlensing experiments, they would increase the 
abundance of heavy metals via the evolution of white d warfs 
and Type I Supernova explosions dCanal et al. 199^ . The 
other problem is that for the mass of the MACHOs to be in 
the range proposed by microlensing observations, the ini- 
tial mass function (IMF) of MACHO progenitors of the 
Galactic halo should be different from those of the disc 
jAdams fc Lauehhn 1995 IChabrier et al. 19961) . otherwise 
we should observe at the tail of mass function (MF) a large 
number of luminous stars and heavy star explosions in the 
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Galactic halos. 

In this study we use the hypothesis of a spatially vary- 
ing MF (instead of uniform Dirac-Delta MF for the halo 
MACHOs) to interpret the LMC microlensing candidates 
(Kerins & Evans 1998, hereafter KE). The physical mo- 
tivation for the hypothesis of spatially varying MFs of 
MA CHOs comes from baryonic cluster formation theo- 
ries jAshman 199d:ICarr lOollDe Paolis et al. 19951) . These 
models predict the spatial variation of MF in the galactic 
halo in such a way that the the inner halo comprises partly 
visible stars, in association with the globular cluster popu- 
lation, while the outer halo comprises mostly low-mass stars 
and brown dwarfs. 

We extend the work of KE by (i) usi ng spatially varying MF 
model in the power-law halo model iAlcock et al. 199al . in- 
cluding the contribution disc, spheroid and LMC for com- 
parison with the latest (5 yr ) LMC microlensing data (Al- 
cock et al. 2000) ; (ii) using a statistical approach applied by 
Green & Jedamzik (2002) and Rahvar (2004) to compare the 
distribution of the duration of the observed events with the 
galactic models; and (iii) performing a likelihood analysis to 
find the best parameters of the inhomogeneous MF model. 
The advantage of using spatially varying MF models is that 
the active mean mass of lenses as the mean mass of observed 
lenses is always larger than mean mass of overall lenses. This 
effect is shown by a Monte-Carlo simulation, and taking it 
into account may resolve the problems with interpreting mi- 
crolensing data. 

This paper is organized as follows. In Section 2 we give a 
brief account of the hypothesis of spatially varying MFs and 
the galactic models used in our analysis. In Section 3 we 
perform a numerical simulation to generate the expected 
distribution of events, taking into account the observational 
efficiency of the MACHO experiment. In Section 4 we com- 
pare the theoretical distribution of the duration of the events 
with the observation. We also perform a likelihood analysis 
to find the best parameters of the MF for compatibility with 
the observed data. The results are discussed in Section 5. 



2 MATTER DISTRIBUTION IN THE 
GALACTIC MODELS 

Spiral galaxies have three components: the halo, the disc 
and the bulge. We can combine these comp onents to build 
various galactic models jAlcock et al. 19961) . In this section 
we give a brief account on the power-law halo and disc mod- 
els which can contribute to the LMC microlensing events. In 
the second part we discuss about MFs of MACHOs and our 
physical motivation for considering spatially varying MFs. 



2.1 Power-law halo mode 

A large set of axisymmetric galactic halo models are the 
"power l aw " models with a matter density distribution 
given by llEvans 1994^ : 



p{R,z) 



47rGg2 

R,\l + 2g^) + R^l - I3q^) + z^[2 - (1 + P)/q^] 

JrJTWT^^T'^W^^ 



where R and z are the coordinates in the cylindrical system, 
Rc is the core radius and q is the flattening parameter which 
is the axial ratio of the concentric equipotentials, the param- 
eter P determines whether the rotational curve asymptoti- 
cally rises, falls or is flat and the parameter Va determines 
the overall depth of the potential well and hence gives the 
typical velocities of objects in the halo. The dispersion ve- 
locity of particles in the halo can be obtained by averaging 
the square of the velocity over the phase space. 
Apart from the Galactic halo, there are other components 
of the Milky Way such as the Galactic disc, spheroid and 
LMC disc that can contribute to the LMC microlensing 
events. The matter distribution of the disc is described by 
double exponentials iBinnev fc Tremaine 198^ and the MF 
of this structure is taken according to the Hubbel Space 
Telescope (HST) observations (Gould., Bahcall & Flynn 
1997). The second component of the Milky Way which 
may also contribute to the microlensing events is the Milky 
Way Spheroid. The spheroid density is given by pspher = 
1.18 X lO"''(r/_Ro)"^ '^M0pc"^, where i?o is the distance of 
the Sun from the center of Galaxy (Guidice., MoUerach & 
Roulet 1994; Alcock et al 1996). We take the dispersion ve- 
locity for this structure to be (t„ = 120km/ s. The mat- 
ter density of LMC disc as the third structure that can 
contribute to the microlensing events is also described as 
a double exponential with the parameters Rd = 1.57kpc, 
h = O.Sfcpc and a — 25km/ s, where Rd is the disc scale- 
length, /!_Js_th«_dlsc_scalehight and a is the dispersion ve- 
locity llGvuk et al. 200dl . 

The combination of galactic substructures as galactic mod- 
els denoted by S, A, B ,C', D, E, F and G. The parameters of 
these models are given in Table. 1. 



2.2 Spatially varying mass function 

The tradition in the interpretation of gravitational mi- 
crolensing data is to use the Dirac-Delta function as the 
simplest MF of Galactic halo MACHOs. Colour-magnitude 
diagram studies of the population of stars in the Galactic 
disc, bulge and other galaxies show that MF behaves like 
a power law function, where the mean mass of the stars 
depends on the density of interstellar medium where the 
stars have been formed. Following this argument, the MF of 
MACHOs in the Galactic halo may also follow a power-law 
function. Fall & Rees (1985) proposed a cooling mechanism 
for the globular cluster formation and on the same basis, 
the hydrogen clouds cooling mechanism can produce a clus- 
ter of brown dwarfs (Ashman 1990). The dependence of the 
mass of stars on the density of the star formation medium 
may causes the heavy MACHOs produced in the dense inner 
regions and the light ones at the diluted areas of the halo 
boarder. 

Here in our study we use MF(r) — 5[M — M{r)] as 
the simplest spatially varying MF, proposed by KE. The 
mass scale M(r) in this model decreases monotonically as 
Mu{^^ Y^'^'^'^'-° , where r is the distance from the center of 
the Gafaxy, Ml and Mu are the mass scales representing 
the lower and upper limits of the mass function and Rhaio 
is the size of Galactic halo contains MACHOs. Considering 
the cold dark matter component for the the halo, the size of 
the Galactic halo may be larger than Rhaio- 
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(1) 


Description 


Medium 


Medium 


Large 


Small 


E6 


Maximal 


Thick 


Thick 


(2) 


/3 







-0.2 


0.2 














(3) 


9 




1 


1 


1 


0.71 


1 


1 


1 


(4) 


Va{km/s) 




200 


200 


180 


200 


90 


150 


180 


(5) 


Rc{kpc) 


5 


5 


5 


5 


5 


20 


25 


20 


(6) 


Ro(kpc) 


8.5 


8.5 


8.5 


8.5 


8.5 


7 


7.9 


7.9 


(7) 


So (M0/pc2) 


50 


50 


50 


50 


50 


100 


80 


80 


(8) 


Rd(kpc) 


3.5 


3.5 


3.5 


3.5 


3.5 


3.5 


3 


3 


(9) 


h{kpc) 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


1 


1 


(10) 


<Tv{km/s) 


31 


31 


31 


31 


31 


31 


49 


49 



Table 1. The parameters of the eight Galactic models: First line is the description of the models in terms of the disk, second line the 
slope of rotation curve (/3 = flat, /3 < rising and /3 > falling), third line the halo flattening (q = 1 represent spherical), fourth line 
(va) the normalization velocity, fifth line Rc halo core, sixth line distance of the sun from the center of galaxy, seventh line the local 
column density of the disk (Eq = 50 for canonical disk. So = 80 for maximal thin disk and So = 40 for thick disk), eighth line disk 
scalelength, the ninth line disk scalehight and tenth line is the adopted one-dimensional velocity dispersion of disk, perpendicular to our 
line of sight. 



3 MICROLENSING EVENTS IN THE 
SPATIALLY VARYING MF 

In this section our aim is to generate microlensing events 
in the spatially varying MF and compare them with the 
observed data. The overall rate of microlensing events owing 
to the contribution of the halo, disk, spheroid and LMC itself 
is given by 

^ = f^lhalo) + —(disk) + — (spheroid) + — {LMC),(2) 
at ' at at at at 

where, / is the fraction of the halo made by the MACHOs. 
The parameter / can be obtained by comparing the ob- 
served optical dept h with that of the expected value from 
the galactic models jAlcock et al. 1995^ ■ The observed opti- 
cal depth is given by Tabs = ^Sii , for 13 microlensing can- 
didates of the MACHO experiment (see Table 3) where E = 
6.12 X 10^ objects- years exposure time. Tabs = 4.43 x 10~* 
iAlcock et al. 200(11 . The observed optical depth is sensitive 
to our estimation of the duration of the events. On the other 
hand, the theoretical optical depth is given by 



'expected 



dT 
HI 



e(t)tdt, 



(3) 



where e{t) is the observational efficiency. Table 2 shows the 
results of comparison between the theoretical and observed 
optical depths. We use the evaluated value of / in each model 
to obtain the distribution of the duration of events. Fig. 
compares the normalized distributions of the duration of 
events for uniform and spatially varying MFs in eight galac- 
tic models. 

The advantage of using a spatially varying MF is that heavy 
MACHOs in the Galactic halo contribute to gravitational 
microlensing events more frequently than do dominant-light 
MACHOs. This effect can be shown by a Monte-Carlo sim- 
ulation. Before explaining the simulation we introduce two 
parameters of the passive and active mean masses of lenses. 
The passive mean mass is defined as the mean mass of the 
overall lenses of the Galactic halo. This mass is independent 
of the gravitational microlensing observation and can be ob- 



tained directly by averaging over the masses of MACHOs 
/ (l>[M, x]M{x)d^xdM 



< M >„ 



J (l>[M, x]d3xdM 
J p(x)d^x 

J M(x)"- 



(4) 



the second equation is obtained by substituting the spatially 
varying MF model. 

In contrast to the passive mean mass, we define the active 
mean mass of lenses as the mean mass of observed microlens- 
ing candidates. It is clear that in the case of a uniform Dirac- 
Delta MF, these two masses are equal, but in the spatially 
varying MF, the active mean mass of lenses is always larger 
than the passive one. 

The algorithm of our simulation for evaluating the active 
mean mass of lenses is (i) selecting the position of lenses 
according to the position distribution function of MACHOs 
along our line of sight; (ii) calculating duration of the events 
(te) and comparing them with the observational efficiency of 
MACHO experiment and (iii) calculating the mean mass of 
selected evens. To select the location of a lens, we imagine 
that we make observations for a given interval of Tabs- The 
probability that a MACHO is located at a distance x = ^ 
from the observer playing the role of a lens, thereby magni- 
fying one of the background stars of the LMC, is 



dT 
dx 



GDs 
M{x)c- 



:x(l - x)vt{x)p{x) 



(5) 



where M{x) is the mass of the MACHO and can be 
substituted by the spatially varying MF model and Vt{x) is 
the transverse velocity of the lens with respect to our line 
of sight. The duration of the events (after picking up the 
location of lenses) is obtained by te — 



vt{x) 



Each time 

at the the Monte-Carlo simulation loop, by comparing the 
duration of the event with the observational efficiency of 
the MACHO experiment, the event is selected or rejected. 
The mass of selected events are used to calculate the mean 
mass of observed MACHOs. Table |5] shows the results of 
our simulation, the passive < Mmi > and active < Mmi > 
mean mass of lenses for different galactic models. As we 
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Figure 1. The normalized distribution of events duration in the 
galactic models of Standard, A, B, C, D, E, F and G, multiplied 
to the observational efficiency of MACHO experiment. The solid 
and dashed lines represent the distributions of events in the uni- 
form and spatially varying Dirac-Delta MFs. The dash-dotted 
line is resulted from the likelihood analysis for the best parame- 
ters of KE model to be compatible with the observed data. The 
duration distributions indicated by solid and dashed-dotted are 
similar while the dashed-line which represents the KE model is 
different than them. 



Figure 2. The expected distributions of the mean (left column) 
and the width (right column) of the duration of events are shown 
for the Standard, A, B and C Galactic models. The cross in- 
dicates the values of < > and h.te from the 13 microlensing 
candidates of MACHO experiment. These distributions arc shown 
for three categories of uniform Dirac Delta MF (solid line), KE 
spatially varying MF (dashed line) and KE model with the pa- 
rameters derived from the likelihood analysis (dash-dotted lines) . 
Most of the galactic models in the spatially varying MF with the 
optimized parameters are compatible with the observed data. 



expected, in all the galactic models the active mean mass is 
always larger than the passive one. This means that in spite 
of the hght abundant brown dwarfs in the Galactic halo, 
lenses with the larger masses produce most microlensing 
events. 



4 COMPARISON OF THE MICROLENSING 
CANDIDATES WITH THE GALACTIC 
MODELS 

In this section our aim is to compare the expected events 
from the spatially varying MF with the microlensing can- 
didates. The next step is to find the best parameters for 
the spatially varying MF model which are compatible with 
the data. Two statistical parameters, the width of the dis- 
tribution of the duration of the events and its mean value 
are used i n our comparison. These parameters a re defined 
as follows JGreen fc Jedamzik 2o6liRahvar 2004^ : 

Ate = Max{te) - Min{te), (6) 
<U> = 1/iVEte. (7) 

Ate and < > for the LMC candidates are 188 and 97 
d, respectively (see Table O. We perform a Monte Carlo 



simulation to generate the mentioned statistical parameters 
from the theoretical distribution of te for comparison with 
the observations. In this simulation we make an ensemble 
of 13 microlensing events where those events are picked up 
from the theoretical distribution of the duration and in each 
set of events, the mean and the width of the duration of the 
events are calculated. The mean of and the width from each 
set is used to generate the distributions. 

With this procedure we obtain the distributions of Ate 
and <te> for three categories of (i) Dirac-Delta MF; (ii) a 
spatially varying MF; and (iii) a spatially varying MF with 
the optimized parameters compatible with the data, result- 
ing from the likelihood analysis. Figs. Inland |H] compare the 
distributions of the observed Ate and <te> with three MF 
models used in eight power-law galactic models. Comparing 
the observed value, indicated by cross in Figs|21and|31 with 
the theoretical distributions of < te > and Ate, shows that 
for Dirac-Delta MF, some of the galactic models such as A, 
C and E are in agreement with the observations while for 
the KE model non of them are compatible with the data. 
To find the spatially varying MF model that is compatible 
with the observations, we perform a likelihood analysis to 
obtain the best upper limit of the MACHO mass and the 
size of the halo in the KE model. The results of analysis in 
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Events{l) Model(2) MF{3) halosize{kpc){i) Ml{5) A4u{6) < M^j > (7) < Af„; > (8) fML(9) 
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Table 2. The first column gives the number of mierolensing events that have been observed during 2 or 5.7 yrs monitoring of LMC 
stars by the MACHO group. The second column indicates the name of eight galactic models as described in Table 1. The third column 
specifies the MF in each model where U indicates the uniform Dirac Delta MF which has been obtained by MACHO group, KE indicates 
the MF proposed by Kerins and Evans (1998) and LA indicates the MF which has been obtained by the likelihood analysis. The fourth 
column shows the size of halo that MACHOs are extended. The fifth column is the lower limit for the mass of MACHOs that are located 
at the edge of halo and the sixth column is the upper limit for the mass of MACHOs that reside at the center of halo. The seventh 
column is the mean mass of the MACHOs in each model, so-call passive mean mass of the lenses and the eighth column is the active 
mean mass of the observed lenses by the experiment. The ninth column shows the halo fraction made by MACHOs in each model. 



Event ; 1 4 5 6 7 8 13 14 15 18 21 23 25 

tE (days) 34.5 83.3 109.8 92 112.6 66.4 222.7 106.5 41.9 75.8 141.5 88.9 85.3 



Table 3. Microlensing candidates observed by the MACHO experiment during 5.7 yrs of observing 11.9 million LMC stars (Alcock et 
al. 2000). First row gives the name of the event according to numbering used by the MACHO group and second row shows the duration 
of event. 



the power-law galactic models are shown in Table|5|with the 
corresponding distribution of the duration of events in Fig. 
1. Figs. |5| and |3| show that the standard model and models 
A, B, C and D using this MF are in good agreement with 
the data. 

In addition to the hypothesis of a spatially varying MF, 
there are other hypothesis, such as self-lensing, that need 
to be confirmed using sufficient statistics of microlensing 
events. Recent microlensing surveys such as Optical Gravi- 
tational Lensing Experiment (OGLE) 3 and SuperMACHO" 
are monitoring LMC stars and will provide more microlens- 
ing candidates over the coming years. To use the results of 
our analysis in the mentioned experiments, we obtained the 
theoretical distribution of events in each model without ap- 
plying any observational efficiency (Fig. |1J. The expected 

3 http: / /bulge.princeton.edu~ogle/| 

* http: //www.ctio.noao.edu~supermacho/ 1 



distribution of events in each experiment can be obtained 
by multiplying the observational efficiencies to these theo- 
retical models. 

It is worth to mention that our statistical analysis is sen- 
sitive to our estimation of the duration of the microlensing 
candidates. The correction with the blending effect can al- 
ter our result. The blending effect makes a source star to be 
brighter than its actual brightness, and the lensing duration 
appears shorter. The duration of a microlensing event can 
be determined from a light curve fit in which the bright- 
ness of the source star is included as a fitting parameter. 
The main problem with this standard method is the de- 
generacy caused by the fitting. High-resolution images by 
the HST have been used to resolve blen ding by random fiel d 
stars in eight LMC microlensing events jAlcock et al. 200 ll) . 
The MACHO group also used another procedure where each 
event is fitted with a light curve that assumes no blending 
and then a correction is applied to the time-scale to account 
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Figure 3. The expected distributions of tlie mean (left column) 
and the width (right column) of the duration of events are shown 
for the D, E, F and G Galactic models. The cross indicates the 
values oi < te > and Ate from the 13 microlensing candidates 
of MACHO experiment These distributions are shown for three 
categories of uniform Dirac Delta MF (solid line), KE spatially 
varying MF (dashed line) and KE model with the parameters 
derived from the likelihood analysis (dash-dotted lines). Most of 
the galactic model in the spatially varying MF with the optimized 
parameters, except model E are compatible with the observed 
data. 



for the fact that blending tends to make the time scales 
appear shorter. This correction was determined from the ef- 
ficiency of a Monte Carlo simulation, and it is a function 
of the time-scale of the measured event. The procedure is 
designed to give the correct average event time-scale, but 
it do es not preser ve the width of the time-scale distribu- 
tion l^enngOflS)- Green & Jedamzik (2002) and Rahvar 
(2004) showed that the width of duration of events derived 
from this method is narrower than the theoretical expecta- 
tions. 



5 CONCLUSION 

In this work we extended the hypothesis of a spatially 
varying MF proposed by KE as a possible solution resolv- 
ing discrepancies between microlensing results and other 
observations. The main point of this is to investigate the 
contradiction where microlensing experiments predict large 
numbers of white dwarfs which have not been observed. 
The advantage of using a spatially varying MF is that we 
can modify our interpretation of microlensing data. We 
showed that in this model, in contrast to the Dirac-Delta 
MF, massive MACHOs contribute in the microlensing 




Standard Model 
Model A 
Model B 
Model C 
Model D 
Model E 
Model F 
Model G 





Figure 4. The duration distribution of events in eight galactic 
models with the spatially varying MF, where the MF parameters 
are optimized to be compatible with the data by the likelihood 
analysis. The expected distribution of events can be obtained by 
multiplying the duration distribution to the observational effi- 
ciency of microlensing experiment. 



events more frequently than the low-mass ones do. To 
quantify our argument we defined two mass scales, the 
active mean mass of MACHOs as the mean mass of lenses 
that can be observed by the gravitational microlensing 
experiment and the passive mean mass of MACHOs as 
the overall mean mass of them. We showed that the active 
mean mass of MACHOs is always larger than the passive 
mean mass, except in the case of a uniform Dirac Delta MF 
where they are equal. 

To test the compatibility of this model with the observed 
microlensing events, we compared the duration distribution 
of the events in this model with the LMC candidates of 
MACHO experiment. We used two statistical parameters 
- the mean and the width of the duration of events - to 
compare the observed data with the theoretical models. 
We showed that amongst power-low halo models some 
of them with a Dirac-Delta MF are compatible with the 
data, while in the case KE model, almost none of them 
are compatible with the data. The best parameters for 
the KE model were obtained with likelihood analysis. In 
the spatially varying MF using the new parameters, some 
Galactic models (such as standard model and models A, B, 
C and D) were compatible with the data. The hypothesis 
of a spatially varying MF of MACHOs may be tested 
by measuring the proper motions of white dwarfs in the 
Galactic halo (Torres ct al. 2002.) . 
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